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Recent thermal conductivity measurements on U Pts single crystals by Lussier et al. indicate 
the existence of a strong b-c anisotropy in the superconducting state. We calculate the thermal 
conductivity in various unconventional candidate states appropriate for the U Pts "B phase" and 
compare with experiment, specifically the i?2ii and Eig states predicted in some Ginzburg- 

Landau analyses of the phase diagram. For the simplest realizations of these states over spherical or 
ellipsoidal Fermi surfaces, the normalized i?2ii conductivity is found, surprisingly, to be completely 
isotropic. We discuss the effects of inelastic scattering and realistic Fermi surface anisotropy, and 
deduce constraints on the symmetry class of the U Pts ground state. 
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The heavy fermion superconductor U Pts is now gen- 
erally thought to be "unconventional" , in the sense that 
the order parameter exhibits a lower symmetry than that 
of the Fermi surface. The principal evidence for this con- 
clusion has been the observation of a nontrivial phase 
diagram for the system in applied magnetic field and 
pressure. Earlier hints of unconventional behavior in- 
cluded transport and thermodynamic properties reported 
to vary as power laws in the temperature for T <t^Tc, as 
well as strong, temperature-dependent anisotropy in su- 
perconducting transport properties. [1] 

Current Ginzburg-Landau (GL) theories of the U Pts 
phase diagram attempt to explain the multiplicity of 
superconducting phases observed by assuming an order 
parameter corresponding to 1) a higher dimensional ir- 
reducible representation of the normal state symmetry 
group, weakly split by a symmetry-breaking field; or 2) 
a mixture of two accidentally nearly degenerate represen- 
tations. In most scenarios, the quartic terms in the GL 
free energy are chosen to reproduce, if possible, the ex- 
pected nature of the low-temperature, low-field U Pts "B- 
phase" , characterized by a line of gap nodes in the hexag- 
onal basal plane, and possibly point nodes along the c- 
axis. The prejudice in favor of this structure arises from 
early tranport experiments. [1,2] Qualitatively speaking, 
these experiments support a picture in which there is, for 
T Tc, a higher density of excited quasiparticles with 
wave vectors in the basal plane. Attempts to measure 
this anisotropy in other experiments have not been uni- 
formly successful, however. For example, measurements 
of the London penetration depth in geometries cho- 
sen to maximize currents along chosen directions found 
AAi(T) ~ with coefficients depending only weakly on 
direction. [3] A /^SR experiment reporting considerable 
anisotropy in the same quantity [4] may in fact have been 
dominated by extraneous effects. [5] Explicit evidence in 
favor of strong gap anisotropy in the superconducting 
state involving measurements in different directions in 
the same sample is in fact rather limited. [1,6] The re- 
cent thermal conductivity measurements of Lussier et al. 
are therefore of considerable interest, both as hard quali- 



tative evidence for anisotropy and for the opportunity to 
make quantitative comparisons with theory. 

The purpose of this paper is to apply the theory of 
heat conduction in a weak-coupling BCS superconductor 
[7] to the data of Lussier et al., [8] assuming unconven- 
tional candidate order parameters of various types cur- 
rently proposed as ground states for the U Pts system. 
While it has been argued that fits to transport coeffi- 
cients at low temperatures are insufficient to fix the or- 
der parameter symmetry, given uncertainties in the mod- 
elling of the scattering amplitude, anisotropy ratios of 
transport quantities should depend much more weakly on 
these details, thereby refiecting more directly the actual 
anisotropy of the superfiuid condensate. We therefore fo- 
cus primarily in what follows on the ratio Kc/ki, between 
the conductivities measured for heat currents directed 
along the c- and b-axes, respectively. 

Normal state. A good deal of important informa- 
tion regarding normal state scattering parameters was 
extracted by Lussier et al. from their measurements 
above the critical temperature Tg = 0.5A'. In this region 
their data could be well fit by the form KN,i~^{T)lT = 
(flj- -|- hiT^), where i = h,c indexes the eigenvalues of 
the thermal conductivity tensor. The authors argue con- 
vincingly that i) the heat conductivity they measure is 
almost entirely electronic in origin, noting that a substan- 
tial phonon contribution would not allow for the observed 
isotropy of the Lorenz number np/T. Furthermore, ii) 
both anisotropy ratios aija^ and &j/&c are empirically 
found to be approximately 2.7, suggesting strongly that 
the anisotropy in the electronic thermal conductivity 
Kjv = CyVpT/3 arises from anisotropy in the Fermi ve- 
locity rather than in the relaxation time r (Here Cy is 
the electronic specific heat and vp the Fermi velocity). 
Since vp ^ m~^, this is consistent with the mass ratio 
mii/nic = 1.7 deduced from upper critical field measure- 
ments. [9] We therefore assume i) that heat conduction 
occurs entirely through electronic transport, and ii) re- 
strict our treatment of collisions to s-wave scattering. 

Impurity scattering is treated within a self-consistent 
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t-matrix approximation [10,11] for the self-energy, which 
in the normal state reduces to an s-wave scattering rate 
r, while for inelastic processes we simply broaden all 
quasiparticle states by a phenomenological inverse life- 



time r- 



\T . This crude model will suffice for 



purpose of fixing the relative importance of inelastic con- 
tributions to superconducting state transport. The coef- 
ficient a may be determined from a fit to the Lussier et 
al. normal state coefficient h to satisfy 0.5a ~ T/T^; 
thus at the critical temperature the inelastic and elastic 
scattering rates are roughly comparable. Thus at least 
for the temperature range close to Tg, inclusion of in- 
elastic processes is essential if a quantitative compari- 
son with experiment is to be attempted. The measured 
residual linear term a in k~^^(T) now allows for a crude 
determination of the magnitude of the impurity scatter- 
ing rate. Using the usual kinetic theory result for the 
thermal conductivity, kjv = Cv(Tc)vpT/3, and using the 
linear specific heat coefficient 7 = 0.42J/A'^ • mole and 
the de Haas-van Alphen measurements of Taillefer and 
Lonzarich, [12], we extract T/Tc ~ 0.2 ± 0.1. This is 
substantially larger than impurity scattering rates which 
were earlier considered typical of nominally pure samples, 
[10,13,14] but previous estimates as to the size of this 
quantity were not been based on reliable normal state 
information. 

Superconducting order parameter. As discussed above, 
there is some evidence for an order parameter corre- 
sponding to the low-temperature "B-phase" of UPts 
whose orbital structure manifests line nodes in the hexag- 
onal basal plane. We therefore focus our attention on 
two such states, the Eig (1,«) state proposed by var- 
ious authors [10,15] and the £"211 (1,0 state with spin 
axis (i II c discussed more recently by Norman [16] and 
Sauls [17]. Although the former is an even-parity ("d- 
wave") spin singlet state, and the latter an odd-parity 
("f-wave") spin triplet state, the spin structure of the 
two order parameters is irrelevant for the calculation of 
the thermal conductivity, and will be neglected here. The 
orbital parts of the order parameters At then have sym- 
metry k;,(kx + iky)", where n = 1(2) for the Eig{E2u) 
state. Both states have a line of nodes in the basal plane 
and point nodes along the c-axis but differ significantly in 
the behavior of the order parameter in the vicinity of the 
point nodes. This difference, usually ignored, can have 
important consequences for the anisotropy of transport 
coefficients, as we show below. In the ID (orbital) rep- 
resentation put forward by Machida et al. [18], the most 
likely states are the A'^u state, which supports quasipar- 
ticle excitations with spectrum identical to the £"211 state 
given above, or the polar state with Aj, a k^,. In what 
follows we will therefore refer to Eig, E'^u or polar, mean- 
ing however the broader classes of states with the orbital 
structures given above. 

The exact form of the order parameter will depend 
on the detailed form of the Fermi velocity ^fv. = 
de{k) I d\i\ip and angle-resolved density of states, as well 



as the distribution of the pair potential weight over the 
various Fermi surface sheets. For the moment we ignore 
all such complications, and assume that the candidate or- 
der parameters in question are represented over a spher- 
ical Fermi surface by At = Ao(T) • sin 20 e"^ for the Eig 

case and At = ^Ao(T) • sin^ 6 cos 6 e^'^ for the E-zu 
case, where Ao(T) = max{|Ai(;|}FS- 

Thermal Conductivity. The thermal conductivity k is 
evaluated using a Kubo formula for the heat-current re- 
sponse similar to the original treatment for an s-wave su- 
perconductor. [7] The generalization to unconventional 
states has been discussed previously by several groups. 
[10,11,19,13] In the limit of vanishing impurity concen- 
trations, identical results were also obtained by Arfi et 
al. [14] using a transport equation method. 

The Kubo formula approach begins with an impurity- 
averaged single-particle propagator which in a Nambu 
matrix notation reads 
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(1) 



where represent the Pauli matrices. Here, we have 
already exploited the assumed particle-hole symmetry 
of the normal state, as well the symmetries of the gap 
functions which lead to vanishing off-diagonal scatter- 
ing self-energy contributions. In this limit, only self- 
energy contibutions to the frequency to, namely Co = 
w — So need to be included. [13] The self-energy So 
due to the elastic impurity scattering is treated in a 
self-consistent T-matrix approximation and is given by 
So = rGo/(c^ — Gq)' where F = nin/^ivNo) is the uni- 
tarity limit scattering rate depending on the concen- 
tration of defects rii, the electron density n, and the 
density of states at the Fermi level Nq- The quantity 
c = cot 60 parameterizes the scattering strength of an 
individual impurity through the s-wave phase shift 60 ■ 
In this work we consider only unitarity limit scattering 
c = since it is clear that weak scattering will lead to 
a weak temperature dependence inconsistent with ex- 
periment for the states in question. [20,21] The inte- 
grated propagator is Gq = (l/2'iTNo)^j^Ti{Tf'g(h,u))}. 
The equation for the self-energies are now solved self- 
consistently together with the gap equation At = 
-3^E.„ Ek' T4k'(l/2)Tr{ri£(k',c.„)} . 

The bare heat current response is given by a convolu- 
tion of the Green's function g with itself at zero exter- 
nal frequency and wave vector weighted with the bare 
heat current vertex wvi^tlf.- [7] Impurity scattering ver- 
tex corrections to current-current correlation functions 
have been shown vanish identically for even parity states 
(Ak = A-k). [13] Even for odd parity states, such cor- 
rections vanish in the unitarity limit considered here. For 
the diagonal thermal conductivity tensor one obtains 
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where u)' and uj" are the real and imaginary parts of uj 
and / is the Fermi function. 

We have evaluated Eqs(2-3) numerically for the two 
candidate states, using impurity scattering rates consis- 
tent with the estimate given above, and compared with 
the data of Lussier et al. [8] If, as we have argued, all nor- 
mal state anisotropy is reflected in the factors of Vp for 
each direction h, c, it is simple to extract the anisotropy 
introduced by pair correlations by plotting the normal- 
ized thermal conductivity (K(T)/T)/(Kff(Tc)/Tc), as in 
Figure 1. We have normalized all quantities to Tc_ since 
the amplitude corresponding to the (1,«) states we con- 
sider first becomes finite there; the comparison with data 
should not be taken seriously in the temperature range 
close to Tc_ , since in a proper theory involving the sym- 
metry breaking field a small admixture of the A-phase 
(1,0) state should be present below the transition. 

We note several interesting aspects of these re- 
sults. Firstly, the £"211 state is found to yield an 
isotropic thermal conductivity within this approxima- 
tion, Kc(T)/ Ki,(T) = 1. This is a consequence of the large 
number of thermally excited quasiparticles in the neigh- 
borhood of the quadratic polar point nodes, which make 
a contribution to the thermal current in both directions 
equal to that of the basal line nodes. The result may be 
shown analytically to hold exactly for all temperatures in 
the superconducting state, for any current-current corre- 
lation function in spherical symmetry. The £"211 state is 
therefore clearly not compatible with the data of Lussier 
et al. in the current approximation (see, however, the 
discussion below). 

The Eig state does exhibit roughly the proper 
anisotropy, as pointed out by Lussier et al., although 
the absolute magnitudes of the calculated thermal con- 
ductivities are too small by roughly 30% at intermedi- 
ate temperatures in this approximation. We note finally 
that at low temperatures T ^ (FTc)^/^ a residual linear 
term k ~ agT should be observable, at least in the b 
direction. [10,11] This "gapless" behavior refiects a cor- 
responding residual density of quasiparticle states in the 
self-consistent treatment of impurity scattering for states 
with line nodes, absent in the theory of Arfi et al. [14] 

Inelastic scattering. The discrepancy between the 
magnitude of the theoretical and experimental results in 
the temperature range near T^, as well as the extracted 
magnitudes of the normal state relaxation components, 
suggests that some treatment of inelastic scattering is 
necessary if a semiquantitative fit is to be attempted. 
In the absence of a complete microscopic description of 
the electron-electron collision processes in this system, 
we simply assume that below the superconducting transi- 
tion, the number of quasiparticles available for scattering 
varies as T/Aq due to the line nodes in the gap. [22] The 
model Tin~^ is therefore a simple interpolation between 



the two limits 



aT2.T/Ao(0) T<Te 
T<T, 



(4) 



We find that our results do not depend significantly on 
the choice of interpolation or the prefactor of T/Aq. 

O 

csj 



h-^ q 



o 
6 




^ r/Tco=0.1 



0.00 0.25 0.50 0.75 1.00 



c- 

FIG. 1. Normalized thermal conductivity for impurity scat- 
tering rate T /Tco = 0.1, no inelastic scattering. Tco is 
transition temperature in absence of impurities. Solid line 
(dashed lines): normalized Hb/T (kc/T) for Eig state. 
Dashed-dotted line: normalized Hb/T and Kc/T for i?2ii state. 
Squares (circles): Kb/T{Kc/T) from reference [8] Insert: vari- 
ation with r. 

In order to get a qualitative picture about the infiu- 
ence of inelastic scattering contributions we now use a 
self-energy So = — F/Go — {i/2)T~^(T), where Go is cal- 
culated as before but using u) = w — So using the modified 
So. In doing so we neglect not only the real part of the 
self-energy but also vertex corrections due to inelastic 
processes. 

In Figure 2, we see that qualitatively, the effect of in- 
elastic scattering is to increase the rise of the normalized 
K/Tjust below Tc as expected. Quantitatively, it is in- 
triguing to note that the inclusion of inelastic processes 
at this level of approximation is sufficient to bring the 
theoretical curves into rough agreement with experiment 
in the intermediate to high temperature range. 

Anisotropy. All the above results depend directly upon 
our assumptions regarding the relaxation time r. The 
ratio of the thermal conductivities Kc/ki, in the normal 
state is consistent with a model in which all anisotropy 
arises through the Fermi velocities. If this is true in the 
superconducting state as well, the anisotropy ratio is only 
very weakly dependent on the relaxation time, and thus 
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constitutes an excellent probe of the intrinsic anisotropy 
of the superconducting condensate. In Figure 3 we plot 
the anisotropy ratio for the data of Lussier et al. and 
compare with our two theoretical candidate states. The 
figure suggests that the effective order parameter struc- 
ture corresponds to a distribution of nodes intermediate 
between the Eig and £"211 model states considered. 
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FIG. 2. Effect of inelastic scattering on thermal conductiv- 
ity. Theory and experimental data as described in Fig. 1, but 
with inelastic self-energy effects included. Dashed-dotted line 
is for i?2ii state, both Kb and Kc. 

The most obvious source of error in the above anal- 
ysis is the neglect of Fermi surface anisotropy, and at 
first glance it would appear that the simplest ellipsoidal 
model accounting for the mass anisotropy mj/mc = 1.7 
observed in U Pts should be sufficient to break the ar- 
tificial degeneracy between the two eigenvalues of k in 
the £"211 state. A straightforward analysis shows, how- 
ever, that this is not the case, and that all ellipsoidal 
anisotropy is refiected in the ratio Vp^/vpi,, which di- 
vides out of the normalized quantity plotted in Fig. 3. 
It is of course clear that inclusion of the complex struc- 
ture of the true U Pts Fermi surface will give rise to some 
anisotropy, but significant effects of this kind appear un- 
likely. Finally, we note that the angular functions con- 
sidered here are merely highly probable members of an 
infinite sequence of basis functions belonging to each ir- 
reducible representation. Inclusion of further such terms 
could also alter the anisotropy. While we can draw no 
definitive conclusions pending analysis of these effects, 
the data do appear to favor a B-phase ground state with 
line nodes in the basal plane and point nodes vanishing 
linearly near the poles. 
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FIG. 3. Thermal conductivtity anisotropy Kc/k6 of various 
states including inelastic scattering. Polar state (Ak = Aokz) 
is included for reference. 
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